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ABSTRACT 


Due  to  increasing  concerns  about  CO2  emissions  and  the  shortage  of  fossil  fuels, 
renewable  energy  has  become  a  major  topic  in  economic  discussions.  One  renewable 
source  is  energy  that  can  be  extracted  from  the  wind.  This  thesis  covers  the  basics  of 
using  a  doubly-fed  induction  generator  (DFIG)  to  convert  the  mechanical  energy  of  the 
wind  into  useful  electrical  power  that  can  be  used  to  supply  electricity  to  any  grid. 
Implementation  and  simulation  results  are  analyzed  in  this  research.  The  design 
implements  digital  four  quadrant  control  of  a  DFIG  with  a  direct  current  (DC)  machine 
serving  as  the  prime  mover.  Digital  control  of  voltage,  current,  and  frequency  in  the  rotor 
windings  is  accomplished  using  a  Voltage  Source  Inverter  (VSI),  while  the  stator  voltage 
and  frequency  is  maintained  by  the  grid.  Simulation  is  accomplished  using  Matlab  and 
Simulink  software.  The  simulations  are  verified  with  lab  hardware. 
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EXECUTIVE  SUMMARY 


Due  to  increasing  concerns  about  CO2  emissions  and  the  finite  supply  of  fossil 
fuels,  renewable  energy  has  become  a  major  topic  across  the  globe.  Wind  energy 
generation  has  attracted  much  interest  in  the  last  few  years.  Large  wind  farms  have  been 
planned  and  installed  in  various  locations  around  the  world.  Many  of  these  wind  farms 
are  based  on  the  Doubly  Fed  Induction  Generator  (DFIG)  technology  with  converter 
ratings  around  30  percent  of  the  generator  ratings  [1],  DFIGs  have  some  advantages  over 
synchronous  and  induction  generators  when  used  in  wind  fanns,  such  as  variable  speed 
operation,  active  and  reactive  power  control,  and  lower  converter  cost  [1]. 

The  objective  of  this  research  is  to  simulate  and  implement  a  wind  energy 
conversion  system  via  a  DFIG.  The  control  strategy  utilizes  a  Field  Programmable  Gate 
Array  (FPGA)  and  a  Voltage  Source  Inverter  (VSI)  to  implement  digital  four  quadrant 
control  of  the  rotor  windings  of  a  DFIG.  Matlab/Simulink  is  the  software  resource  used 
to  design,  simulate,  and  predict  the  system  behavior  before  actual  implementation.  A 
Direct  Current  (DC)  machine  is  used  as  mechanical  load  representing  the  wind  energy 
source  during  the  implementation. 


DFIG 


Figure  1.  Simplified  Circuit  of  DFIG  Control  for  Wind  Energy  Conversion. 


xv 


As  seen  in  Figure  1,  the  voltage  source  inverter  converts  a  3-phase  AC  source  to  a 
controlled  DC  bus  voltage.  The  DC  voltage  is  then  inverted  into  a  3 -phase  wavefonn  via 
a  digitally  controlled  VSI.  The  digital  control  for  the  Insulated  Gate  Bipolar  Transistor 
(IGBT)  network  is  implemented  via  a  FPGA.  The  FPGA  interface  samples  the  voltage 
and  current  of  a  DFIG  and  speed  of  the  DC  drive  through  an  analog  to  digital  (A/D) 
converter.  An  algorithm  then  calculates  the  phase  and  amplitude  of  the  voltage  and 
current  needed  at  the  terminals  of  the  rotor  circuit.  By  controlling  the  phase  and 
amplitude  of  the  currents  in  the  rotor  windings  various  modes  of  operation  and  control  of 
the  DFIG  can  be  accomplished.  Real  power  flow  and  reactive  power  flow  can  be 
controlled,  as  shown  in  Figure  2,  where  the  variables  Vs  and  Is  represent  the  stator  voltage 
and  stator  current. 
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Figure  2.  Power  Factor  Control  of  the  stator  circuit  is  achieved  by  changing  the  rotor 

currents. 
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A  DC  machine  is  used  as  the  prime  mover  to  simulate  the  wind.  By  manually  controlling 
the  torque  produced  by  the  DC  machine,  variable  wind  load  conditions  can  be 
represented.  As  input  torque  is  increased,  the  speed  of  the  rotor  changes  and  the  position 
is  fed  into  the  algorithm  of  the  FPGA.  This  corresponding  position  is  then  used  to  define 
the  current  in  the  rotor  windings.  The  increase  in  input  torque  increases  the  output  power 
produced  by  the  system.  See  variable  torque  operation  in  Figure  3  where  the  stator 
current  increases  when  the  input  torque  increases.  The  variables  Vs,  Is,  and  Ps  represent 
the  stator  voltage,  stator  current,  and  average  instantaneous  output  power  for  one  phase 
multiplied  by  3  from  the  stator. 
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Figure  3.  Step  Change  in  Torque  after  steady  state  is  reached  (stator  current,  stator  voltage 

and  average  instantaneous  stator  power). 
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The  rotor  current  can  also  be  used  to  control  output  power  from  the  generator  during 
steady  state  or  transient  conditions.  An  increase  in  the  rotor  current  during  steady  state 
condition  will  cause  a  change  in  the  stator  current,  which  results  in  a  change  in  the  output 
power.  Figure  4  shows  a  step  change  in  rotor  current  and  how  it  affects  the  output  power 
for  a  given  load  on  the  shaft.  The  output  power  has  slightly  decreased  because  the 
increase  in  rotor  current  has  caused  the  power  factor  to  decrease.  The  variables  I,  and  Ps 
represent  the  rotor  current  and  average  instantaneous  output  power  of  one  phase 
multiplied  by  3  from  the  stator. 
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Figure  4.  Output  Power  vs.  Step  Change  Rotor  Current  after  steady  state  is  reached. 

The  goal  of  this  research  was  to  design  an  electrical  interface  for  a  wind  energy 
conversion  system.  This  thesis  covers  the  simulation  and  hardware  construction  of  a  wind 
energy  conversion  system.  The  simulation  results  will  be  used  to  validate  the  actual 
implementation.  Future  research  will  investigate  control  schemes  to  improve  on 

efficiency,  effectiveness,  and  possibly  grid  fault  analysis. 
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I.  INTRODUCTION 


A.  PURPOSE 

The  mission  at  the  Naval  Postgraduate  School  (NPS)  is  to  provide  students  with 
the  highest  quality  and  most  defense-relevant  graduate  education  available  in  computer 
and  electrical  engineering.  The  United  States  President’s  Agenda  on  Energy  and  the 
Environment  encourages  the  development  of  energy  sources  that  will  increase  national 
security  by  decreasing  our  dependence  on  foreign  oil,  ensure  that  10  percent  of  our 
electricity  comes  from  renewable  sources  by  2012,  and  25  percent  by  2025  [2],  Wind 
energy  is  a  very  promising  renewable  energy  source  of  the  future  and  is  part  of  the 
president’s  overall  strategy  to  develop  a  diverse  portfolio  of  domestic  energy  supplies. 
This  thesis  supports  the  NPS  mission  and  one  topic  discussed  in  the  President’s  Agenda 
on  Energy  and  the  Environment  by  using  a  wind  energy  conversion  system  as  a  means  of 
increasing  national  security  by  harnessing  one  of  nature’s  renewable  energy  sources. 

B.  RESEARCH  OBJECTIVE 

The  goal  of  this  thesis  was  to  build  a  computer  simulation  of  a  Wind  Energy 
Conversion  System,  based  on  mathematical  equations,  through  the  use  of  Mathworks’ 
Simulink  software.  Once  the  simulation  is  complete,  a  fully  functioning  model  of  a  Wind 
Energy  Conversion  System  is  built,  so  that  various  tests  on  the  simulation  and  the 
physical  system  can  be  performed  to  compare  the  simulation  results  with  the  physical 
model.  Based  on  the  successful  results  of  the  system,  we  lay  the  groundwork  for  future 
research  in  the  NPS  power  labs  in  the  field  of  renewable  energy. 

C.  APPROACH 

The  first  step  is  to  generate  a  computer  simulation  of  the  system,  shown  in  Figure 
1.  There  are  numerous  ways  of  simulating  a  DFIG,  but  the  one  used  for  this  thesis  is 
covered  in  [3].  In  particular,  a  computer  model  of  a  DFIG  based  on  mathematical 
equations  is  simulated  in  the  rotor  reference  frame.  A  current  control  algorithm  is 

designed  in  Simulink  along  with  the  Xilinx  toolbox  that  controls  the  rotor  currents  to 

1 


produce  the  required  voltage  and  frequency  needed  on  the  power  grid.  The  wind  is 
simulated  via  a  torque  input  block.  This  model  allows  for  representation  of  various  modes 
of  operation,  which  will  be  useful  for  future  research,  but  for  this  thesis  we  concentrate 
on  variable  load  conditions  on  the  rotor,  which  simulate  changes  in  the  wind  speed.  With 
a  working  simulation  and  the  ability  to  control  rotor  currents,  various  modes  of  operation 
and  control  are  analyzed  through  the  Simulink  software. 

The  physical  system  consists  of  a  DFIG,  VSI,  Student  Design  Center  (SDC)  [4], 
encoder,  and  a  DC  drive  used  to  simulate  changes  in  wind  speed.  The  DFIG  is  connected 
to  the  grid  via  the  stator  windings  and  to  the  VSI  via  the  rotor  windings.  The  SDC  is  the 
mechanism  by  which  the  control  algorithm  is  implemented.  The  SDC  measures  the 
applicable  system  parameters  and  outputs  the  appropriate  control  signals  to  the  VSI  to 
control  the  rotor  currents  in  a  fashion  that  allows  for  a  specific  mode  of  operation.  In  this 
research,  variable  wind  speed  control,  and  active  and  reactive  power  control  are 
demonstrated. 

D.  THESIS  ORGANIZATION 

This  thesis  is  organized  as  follows.  Chapter  II  covers  the  theory  of  induction 
machines.  The  equations  that  represent  voltage,  flux  linkage,  and  reference  frame 
transfonnations  are  discussed.  The  testing  procedure  used  to  measure  and  calculate  the 
machine  parameters  used  for  simulation  are  also  discussed.  Lastly,  Chapter  II  explains 
how  the  voltage  and  flux  equations  and  machine  parameters  are  used  to  model  the 
system.  Chapter  III  is  dedicated  to  the  control  algorithm.  It  explains  how  an  algorithm  is 
formulated  based  on  the  equations  derived  in  Chapter  II  and  are  manipulated  to  control 
the  simulation  and  the  physical  system.  Chapter  IV  displays  the  results  of  the  working 
system  and  compares  the  actual  results  to  the  simulation  results.  Chapter  V  discusses  the 
conclusions  and  the  continuation  research  objectives  for  this  thesis.  Appendix  A,  B,  C 
and  D  contain  Matlab  files,  data  sheets  for  equipment  used,  a  full  schematic  layout  of  all 
Simulink  systems,  and  transformation  derivations  used  in  this  thesis. 
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II.  DFIG  SIMULATION 


A.  MATHEMATICAL  REPRESENTATION  OF  DFIG 

In  order  to  simulate  the  Wind  Energy  Conversion  System  a  proven  method  is 
needed  to  represent  the  characteristics  of  a  DFIG.  The  equations  used  to  characterize  the 
DFIG  used  in  this  research  are  obtained  from  [3].  The  rest  of  this  section  is  devoted  to  the 
mathematical  equations  that  describe  a  DFIG. 

The  voltage  equations  to  represent  a  3 -phase  induction  machine  can  be  expressed 
as  follows, 

[Vabcs]  =  [rs][iabcs]  + P^abcs] 

[Vabcr]  =  [r,][labcr]  + p[Kbcr] 

where  v,  r,  i,  and  k  respectively  refer  to  the  voltage,  resistance,  current,  and  flux  linkage 
of  the  phase  windings.  The  subscripts  a,  b,  and  c  refer  to  their  phase  component.  The 
subscripts  r  and  s  refer  to  the  stator  and  rotor  windings.  The  term  p  represents  the 
derivative  .  The  flux  linkages  shown  in  (1)  for  a  linear  magnetically  coupled  circuit 
are  expressed  as 
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where  L  refers  to  the  inductance  in  the  respective  winding.  The  subscript  sr  represents  the 
mutual  coupling  between  the  stator  and  rotor  windings.  Finally,  to  complete  equations 
that  describe  an  induction  machine  the  winding  inductances  are  defined  as 
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where  the  variables  Lis  and  Lms  are  the  leakage  and  magnetizing  inductances  of  the  stator 
windings,  Lir  and  Lmr  are  for  the  rotor  windings.  The  variable  Lsr  is  the  amplitude  of  the 
mutual  inductance  between  the  stator  and  rotor  windings  and  9r  is  the  angular  position  of 
the  rotor. 

As  a  standard  practice  and  for  convenience,  the  rotor  variables  are  referred  to  the 
stator  windings  by  the  appropriate  turn  ratio  N. 


4 

4 


4 


^ 'abcr  ^ abcr 


N 

2 '  —  s  2 

abcr  abcr 


r  = 


r  4  v 
v4y 


r. 


(6) 


(7) 


(8) 


(9) 


The  magnetizing  and  mutual  inductances  are  associated  with  the  same  magnetic  flux  path 
and  related  by  the  following  equations. 
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By  defining  L'sr=^-Lsr and  L'mr  = 
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L  the  rotor  inductance  variables  can  be 


referred  to  the  stator  windings  as 
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Substituting  (13  and  14)  into  (1  and  2),  the  flux  linkage  and  voltage  variables  referred  to 
the  stator  windings,  can  now  be  expressed  as 
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(16) 


The  voltages  and  flux  linkages  in  (15  and  16)  describe  the  behavior  of  a  DFIG 
and  have  components  that  vary  with  time  and  rotor  position.  A  change  of  variables  is 
used  to  reduce  the  complexity  of  these  differential  equations.  A  change  of  variables  that 
formulates  a  transfonnation  of  balanced  3-phase  equations  to  a  new  reference  frame  may 
be  expressed  as 
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where  the  matrix  Ks  for  stator  variables  and  Kr  for  rotor  variables  is  expressed  as 
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The  variable  6  in  (18  and  19)  is  the  position  of  the  reference  frame  for  the  fabc  variables 
and P  =  6-6r.  The  x  in  (17)  represents  the  particular  frame  of  reference  that  the 
variables  are  transformed  while  /  represents  the  voltage,  current,  or  flux  linkages  in  the 
stator  or  rotor  windings.  In  the  remaining  equations  in  this  chapter  the  absence  of  a 
variable  in  the  place  holder  of  x  represents  an  arbitrary  frame  of  reference. 
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To  transform  between  reference  frames  use  (20)  where  y  represents  the  new  frame  of 
reference.  Utilizing  a  little  trigonometry,  and  performing  the  transformation  in  (17)  on 
each  component  in  the  matrices  on  the  right  hand  side  of  (16),  the  voltage  equations  can 

now  be  expressed  in  the  arbitrary  reference  frame  as 
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v  =  rJqs + b)K + p\s 
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V'dr  =  r'j'ds  ~  ( ®  -  ®,  )  Ks  +  A 
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where  <y(.  is  the  rotational  speed  of  the  rotor  and  co  is  the  speed  of  the  reference  frame 

that  the  variables  are  transformed  in  radians  per  second.  The  next  step  is  to  express  the 
inductance  component  of  the  flux  linkages  in  (21)  as  a  reactance  instead  of  as  an 
inductance.  This  is  convenient  because  in  the  next  section  (Induction  Motor  Test)  the 
actual  machine  parameters  are  measured  and  calculated  in  ohms.  The  inductance  is 
converted  to  reactance  at  a  base  frequency  ( coh )  of  377  radians  per  second  and  (21)  is 
rewritten  as 
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The  flux  linkage  component  of  (21)  now  has  units  of  flux  linkage  per  second  and 
represented  by  the  symbol^ .  For  ease  of  understanding  the  flux  linkages  in  (21)  and  flux 
linkage  per  second  in  (22)  are  compared 
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Equations  (21  and  22)  are  derived  in  this  section  and  mathematically  represent  the 
behavior  of  a  DFIG.  The  next  section  discusses  how  to  measure  and  calculate  the 
parameters  of  the  actual  machine  used  in  this  thesis. 

B.  INDUCTION  MOTOR  TEST 

In  order  to  efficiently  model  the  system  in  Matlab  and  Simulink,  the  physical 
machines  parameters  are  used  for  the  simulation.  The  DFIG,  model  823 1 ,  used  for  this 
research  is  part  of  the  Lab-Volt  Electro-Mechanical  System  (EMS)  and  consists  of  a  4- 
pole,  0.25  horsepower,  120  volt,  60  hertz,  wave  wound  rotor  machine.  The  stator  consists 
of  a  set  of  3-phase  wye-connected  windings  with  turns  Ns  and  resistance  Rs.  The  rotor 
consists  of  a  set  of  3 -phase  wye-connected  windings  with  turns  Nr  and  resistance  Rr . 
The  induction  motor  test  performed  is  carried  out  as  in  [5], 
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The  steady  state  operating  characteristics  of  the  DFIG  was  analyzed  using  a  per- 
phase  equivalent  circuit  shown  in  Figure  5.  The  DFIG  was  converted  to  an  Induction 
machine  for  the  parameter  test  by  short  circuiting  the  rotor  winding  tenninals  in  a  wye- 
configuration.  The  variables  Vs ,  Rs  and  Xls  represent  the  voltage,  resistance,  and  leakage 

reactance  of  the  stator.  The  variables  Vr\R',  and  X\r  represent  the  voltage,  resistance  and 

leakage  reactance  of  the  rotor  referred  to  the  stator.  Xms  represents  the  stator  magnetizing 
reactance  and  s  is  the  slip. 


Rs 


Xls 


Rr' 


xir 


R[ 


1  -s 

s 


Figure  5.  Induction  Motor  Equivalent  Circuit  (From  [5]). 

1.  DC  Test 

The  first  parameter  test  ran  was  to  calculate  the  stator  winding  resistance  for  each 
phase.  This  test  consisted  of  applying  a  known  DC  voltage  to  the  stator  tenninals  and 
measuring  the  current  flow  through  the  windings.  The  DC  test  on  the  stator  windings  is 
useful  because  it  does  not  induce  a  voltage  in  the  stator  or  rotor  windings  and  eliminates 
the  need  to  measure  or  calculate  any  reactance  caused  by  induced  voltages.  The  only 
circuit  parameter  limiting  current  flow  is  The  measured  voltage  and  current  values  are 

inserted  into  (25)  and  the  stator  resistance  Rs  is  calculated.  This  test  is  performed  on  each 
phase  of  the  stator  windings  and  the  three  values  are  averaged  together. 

Rs= ^  (25) 
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2. 


No  Load  Test 


A  no-load  test  is  performed  to  measure  the  rotational  losses  and  other  equivalent 
circuit  parameters  needed  in  the  blocked  rotor  test.  In  this  test,  rated  voltage  and 
frequency  is  applied  to  the  stator  while  the  machine  is  running  at  no  load.  The  input 
power,  phase  voltage,  and  phase  current  are  measured  for  each  phase  and  averaged.  At  no 
load,  the  slip  is  very  small  and  the  rotor  branch  impedance  is  very  large.  Therefore,  all 
the  real  power  is  dissipated  in  the  stator  resistance.  The  machine  at  no  load  is  very 
inductive,  so  the  input  impedance  is  expressed  as 


V 


a,nL 


=  V^2+(X™+X/J 


(26) 


From  this  equation,  we  can  solve  for  the  quantity  Xm  +  Xls ,  which  is  used  in  the  blocked 
rotor  test. 


3.  Blocked  Rotor  Test 

The  blocked  rotor  test  is  perfonned  to  calculate  the  remaining  circuit  parameters. 
In  this  test,  the  rotor  of  the  induction  machine  is  blocked  and  a  reduced  voltage  is  applied 
to  the  tenninals  so  that  rated  current  flows  through  the  stator  windings  of  the  machine. 
The  input  power,  voltage,  and  current  are  measured  and  averaged.  When  the  rotor  is 
blocked,  the  slip  is  equal  to  one.  As  a  result,  the  rotor  branch  impedance  is  very  small  and 
the  magnetizing  impedance  is  ignored.  Real  power  is  only  consumed  by  the  stator  and 
rotor  resistances.  The  sum  of  the  stator  and  rotor  resistances  can  be  solved  as 


R+R' 


(27) 
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The  input  impedance  is  expressed  as 


V. 
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(28) 


where  it  is  assumed  that  Xls  =  Xlr .  Substituting  the  result  of  (27)  into  (28)  allows  us  to 
solve  for  the  numerical  quantities  of  Xls  and  Xlr .  From  this  point  Xms  can  be  calculated 

by  substituting  the  stator  leakage  reactance  value  into  equation  (26).  Table  1  below  is  a 
summary  of  parameters  calculated  from  the  DFIG  used  for  this  thesis. 


Parameters  (Cl) 

Stator  Resistance  ( r  ) 

12 

Stator  Leakage  Reactance  (A^) 

9.1 

Rotor  Resistance  (/) 

15 

Rotor  Reactance  (X'lrJ 

9.1 

Magnetizing  Reactance  (XnB) 

126 

Table  1.  DFIG  Measured  Values  at  60  FIz. 

C.  DFIG  COMPUTER  MODEL 

The  equations  (22  and  24)  are  the  basis  for  simulating  the  DFIG.  A  block  diagram 
of  the  simulation  setup  is  shown  in  Figure  6.  The  model  simulates  the  DFIG  in  the  rotor 
reference  frame.  Although  the  equations  that  represent  a  DFIG  have  been  defined,  they 
still  need  to  put  in  a  form  that  is  useful  for  computer  simulation.  We  also  need  to  define 
the  torque,  power,  and  rotor  speed  equations,  which  are  needed  for  simulating  the 
machine. 
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Figure  6.  DFIG  Computer  Model  in  Block  Form  (After  [3]). 

The  first  step  is  to  solve  the  flux  linkage  per  second  equations  of  (22)  for  current. 
We  do  this  by  defining  four  new  variables  in  terms  of  quantities  calculated  in  the 
induction  motor  test 
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The  currents  are  now  defined  as 
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Now  that  the  current  equations  have  been  defined,  the  flux  linkage  per  second  equations 
in  (22)  are  rewritten  as  integral  equations  defined  in  terms  of  the  four  new  variables  from 
(29-31).  The  final  form  of  the  flux  linkage  equations  that  are  used  to  simulate  the  DFIG 
shown  in  Figure  6  are  represented  as 
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The  rotor  position,  which  is  needed  to  transform  the  equations  to  the  respective 
frame  of  reference,  is 

(T,-Tl)P 


CO,.  = 


2  J 


dt 


(44) 


where  Te  and  7)  ,  represent  the  electrical  torque  generated  in  the  machine  and  input  torque 

on  the  shaft  of  the  rotor.  The  variables  P  and  J  are  the  number  of  machine  poles  and  the 
inertia  of  rotor.  Once  the  speed  of  the  rotor  is  obtained,  the  position  is  calculated  as 


6r  =  J  cor  dt. 

The  electrical  torque  and  the  power  equations  are  defined  as 


(45) 
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D.  SIMULINK  IMPLEMENTATION  OF  DFIG 

The  DFIG  from  Figure  6  is  implemented  with  Simulink  Software.  The  following 
Figures  7  and  8  show  how  (38-46)  are  implemented  using  the  Simulink  software  to 
represent  the  DFIG. 
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Figure  7.  DFIG  Simulink  Simulation  (From  [6]). 

Figure  7  shows  the  blocks  where  the  flux  linkage  per  second  equations  are 
implemented,  and  how  they  are  used  to  compute  the  rotor  and  stator  currents.  The 
respective  equation  number  is  annotated  on  each  block  along  with  the  associated  input 
and  output.  Figure  8  is  a  drill  down  into  the  block  that  represents  (38)  the  stator  flux 
linkage  per  second  on  the  q-axis.  As  seen  in  Figure  7,  there  are  four  inputs  and  one 
output,  which  also  serve  as  an  input  to  the  same  block.  Using  the  basic  mathematical 
operators  in  the  Simulink  library  (38)  is  represented  as  shown  in  Figure  8.  The  remaining 
blocks  for  each  equation  that  represent  the  DFIG  are  built  the  same  way  and  can  be  found 
in  Appendix  C. 
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Figure  8.  Equation  38. 

In  conclusion,  this  chapter  described  the  methods  in  which  to  simulate  the  DFIG 
used  in  this  thesis.  The  next  chapter  will  discuss  the  methodology  used  to  control  the 
Wind  Energy  Conversion  System. 


16 


III.  CONTROLLER 


A.  CONTROLLER  TOPOLOGY 

As  stated  earlier,  one  of  the  components  of  the  wind  energy  conversion  system  is 
the  rectifier  and  VSI.  The  VSI  and  rectifier  used  for  this  research  are  a  combined  package 
manufactured  by  SEMIKRON  and  called  a  SEMISTACK.  The  SEMISTACK  consists  of 
a  three  phase  rectifier  and  three  phase  inverter.  The  inverters  inside  the  SEMISTACK  are 
made  of  SKM  50  GB  123D  IGBTs  that  are  controlled  by  SEMIKRON  SKHI  22  gate 
drivers.  The  VSI  is  used  to  control  the  magnitude,  frequency,  and  phase  angle  of  the 
currents  in  the  rotor.  The  rectifier  produces  a  constant  DC  voltage  to  the  VSIs.  The 
control  for  each  IGBT  in  the  SEMISTACK  is  implemented  through  the  SDC.  The  SDC 
utilizes  a  Xilinx  Field  Programmable  Gate  Array  (FPGA)  chip  and  a  control  board  that 
has  various  inputs  and  outputs  that  interface  with  the  SEMISTACK  and  DFIG.  One  of 
the  inputs  is  an  8-channel,  12-bit  A/D  converter  used  to  measure  the  voltage  and  current 
signals  from  the  DFIG.  Another  input  samples  the  rotor  speed.  The  pre-programmed 
FPGA  chip,  based  on  the  sampled  inputs,  outputs  the  desired  gate  signals  to  the  control 
board.  The  control  board  has  BNC  (Bayonette  Neil-Concehnan)  connecters  that  connect 
to  the  SEMISTACK’s  gate  drivers.  The  gate  drive  signals  activate  the  IGBTs  in  a  fashion 
that  will  produce  the  desired  3-phase  current  in  the  rotor  windings.  The  SDC  has  many 
parts  and  functions,  which  make  it  useful  as  a  design  and  analysis  tool.  The  procedure  for 
using  the  SDC,  architecture  of  internal  components,  and  the  benefits  of  using  a  FPGA  for 
digital  control  of  power  systems  are  discussed  in  [4]. 

Figure  9  is  a  Simulink  block  diagram  of  the  Wind  Energy  Conversion  System. 
The  diagram  shows  where  the  controller  samples  each  input  and  the  corresponding 
outputs  that  go  to  the  DFIG.  The  controller  is  physically  located  inside  the  SDC  for  the 
hardware  design.  It  is  important  to  note  that  the  same  controller  logic  used  for  simulation 
is  used  in  the  hardware  design;  therefore,  the  input/output  blocks  shown  in  the  simulation 
block  diagram  symbolize  actual  connections  to  the  SDC  from  the  grid  and  the  DFIG. 


17 


Figure  9.  Simulation  of  DFIG  Control  for  Wind  Energy  Conversion. 


Figure  10  is  a  drill  down  of  the  Rotor  Controller  Block.  The  same  logic  blocks 
used  in  the  simulation  to  control  the  DFIG  model  are  used  in  the  hardware  configuration 
to  control  the  DFIG  with  one  exception.  An  additional  set  of  logic  blocks  is  added  to 
implement  Space  Vector  Modulation  (SVM)  which  is  one  of  many  techniques  that  can  be 
implemented  on  a  VSI  to  create  sinusoidal  3-phase  waveforms  and  will  be  discussed 
more  in  Section  B.7.  For  now,  it  is  important  to  know  that  inside  the  rotor  controller  there 
is  a  separate  SVM  block  that  is  not  needed  for  simulation.  Figure  10  shows  how  the  rotor 
controller  is  used  for  the  simulation  and  hardware  design.  The  next  section  explains  how 
each  of  the  components  in  Figure  10  is  used  to  control  the  rotor  current. 
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Figure  10.  Rotor  Controller  Block  Diagram. 


B.  CONTROLLER  COMPONENTS 
1.  Voltage  Source  Simulation 

The  simulation  required  a  voltage  block  to  simulate  the  voltages  on  the  grid.  The 
design  of  the  source  is  shown  in  Figure  1 1 .  The  voltage  block  models  the  rated 
parameters  of  the  DFIG  used  for  this  thesis.  The  voltage  block  produces  three  sinusoidal 

waveforms  at  60  hertz  with  peak  amplitude  of  120*V2and  displaced  equally  120 
degrees  apart.  The  controller  setup  requires  only  two  voltage  inputs  so  the  line  to  line 
values  are  extracted  as  shown. 


Figure  11.  Voltage  Source  Simulation. 
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2. 


Transformation  of  Stator  Voltage 


As  seen  in  Figure  10,  the  stator  voltages  are  fed  into  the  rotor  controller  where 
they  are  transformed  to  the  rotor  reference  frame.  The  transformation  [  fqd0s  1  =  Ks  [  fabcs  ] 
in  (17)  is  slightly  modified  to  transfonn  the  line-to-line  voltages  vabs  and  vbcs  to  the  rotor 

reference  frame  (See  Appendix  D).  Figure  12  demonstrates  how  the  transformation  is 
performed  with  the  Xilinx  library  set.  The  constant  85  in  Figure  12  corresponds  to  a  30 
degree  phase  shift  from  the  transformation  in  Appendix  D.  The  sin  and  cosine  blocks 
from  the  Xilinx  library  use  a  10  bit  lockup  table  to  represent  angles  from  0-360  degrees. 
The  angle  6er  is  obtained  by  taking  the  angle  from  arctangent  of  the  transformed  voltages 

and  represents  the  electrical  position  of  the  currents  in  the  rotor.  This  angle  is  used  to 
transform  the  rotor  currents  to  the  synchronously  rotating  reference  frame.  The  voltage  d- 
axis  component  is  inverted  before  taking  the  arctangent  to  account  for  the  d-axis  being 
1 80  degrees  out  of  phase  from  the  y-axis. 


Figure  12.  Stator  Voltage  Transformation  to  Rotor  Reference  Frame. 

3.  Transformation  of  Rotor  Currents 

The  rotor  currents iabc ,  shown  in  Figure  13,  are  also  fed  into  the  controller  where 
they  are  transformed  using  (17)  to  the  variables/^,,.  This  notation  represents  a 
transformation  of  the  rotor  currents  to  the  synchronously  rotating  reference  frame.  This 
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transformation  produces  a  constant  value,  if  the  currents  are  in  steady  state  condition, 
which  is  useful  in  the  next  block.  The  constant  341.33  represents  a  120  degree  phase  shift 
used  in  (17-19). 


341.33203125 


z  (ab) 


z  Cab) 


z  (ab) 


z  (ab)  -I 


<¥ 


F igure  13.  Rotor  Current  T ransformations . 

4.  Proportional  /Integral  (PI)  Gain 

The  transformation  of  balanced  3-phase  variables  to  the  qdO  axis  will  only 
produce  components  on  the  qd  axis.  Additionally,  if  the  transformation  is  to  the 
synchronously  rotating  reference  frame,  the  q  component  will  have  constant  amplitude 
equal  to  the  amplitude  of  the  phase  variables  and  the  d  component  will  have  amplitude  of 
zero.  This  property  is  useful  in  the  PI  control  algorithm  for  driving  the  steady  state  error 
to  zero.  Figure  14  contains  the  logic  used  for  a  proportional  and  integral  gain  controller 
designed  with  the  Xilinx  library  block  set.  By  transforming  the  rotor  currents  to  the 
synchronously  rotating  reference  frame  two  constants  are  produced  as  discussed  above. 
Setting  the  reference  iq  to  the  desired  output  value  and  the  reference  U  to  zero  allows  for 
full  control  of  the  output.  The  Integral  part  of  the  PI  controller  is  built  with  an 
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accumulator  so  a  reset  is  added  to  set  the  accumulator  to  zero  when  desired.  The  currents 
in  the  rotor  windings  are  controlled  by  adjusting  the  voltage  at  the  rotor  terminals; 
therefore,  the  output  from  the  current  PI  controller  is  a  voltage. 


Figure  14.  PI  Current  Controller. 

5.  Transformation  of  Rotor  Voltages  to  Rotor  Reference  Frame 

The  rotor  voltages  veq  and  Vc/  are  transformed  from  the  synchronously  rotating 

reference  frame  to  the  rotor  reference  frame  as  shown  in  Figure  15.  This  transformation 
utilizes  equation  (20),  which  is  useful  for  transforming  variables  between  frames  of 
reference.  The  input  to  the  VSI  in  Figure  10  is  normalized  to  the  value  of  the  DC  voltage 
produced  by  the  rectifier  on  the  SEMISTACK.  This  normalization  takes  place  in  the  last 
logic  block  of  Figure  15. 
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Figure  15.  Transformation  Between  Reference  Frames. 

6.  Rotor  Encoder 

The  encoder  used  for  this  thesis  is  an  MES20  (Type  C)  incremental  shaft  encoder. 
The  encoder  provides  a  digital  input  to  the  SDC.  The  control  circuitry  inside  the  SDC 
interprets  the  digital  signals  from  the  encoder  and  produces  both  rotor  speed  and  position. 
The  encoder  software  used  in  this  thesis  is  discussed  in  detail  in  [7].  One  modification, 
shown  in  Figure  16,  was  added  to  the  design  in  [7],  which  was  used  to  adjust  rotor 
position  and  speed  for  a  4-pole  machine.  See  Appendix  C  for  complete  encoder  layout. 


From  Encoder 


►  e. 


Figure  16.  Rotor  Encoder  Addition. 
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7.  Space  Vector  Modulation 


One  well-known  way  of  creating  3-phase  waveforms  from  a  VSI  is  via  SVM. 
SVM  is  a  fonn  of  Pulse  Width  Modulation  (PWM)  in  which  an  algorithm  involving 
space  vectors  is  used  to  control  the  on  and  off  times  of  pulsed  signals.  The  pulsed  signals 
then  drive  the  input  signals  to  a  VSI  allowing  the  user  to  create  any  magnitude  and 
frequency  of  signal  desired.  The  SVM  approach  utilized  in  this  controller  is  obtained 
from  [7]. 


Figure  17.  VSI  and  SVM  Hexagon  (From  [7]). 

The  VSI  produces  a  three  phase  output  voltage  that  is  controlled  by  turning  on 
and  off  the  six  IGBTs  on  the  VSI  diagram  in  Figure  17.  The  on  and  off  switching  states 
are  based  on  the  SVM  Hexagon  also  shown  in  Figure  17.  The  qd  axis  is  overlaid  on  the 

SVM  Hexagon  and  represents  the  axis  of  the  variables  Vs  and  vsdr,  which  are  shown 

entering  the  SVM  block  in  Figure  10.  The  magnitude  of  the  vectors  vqr  and  vsdr  is  equal  to 

the  vector  V*  and  0  is  equal  to  the  arctangent  of  the  two  vectors. 

There  are  six  sectors  on  the  hexagon  and  eight  possible  on  and  off  states  that  are 
v*  V  V 

used  to  produce  the  vectors  v  .  The  vectors  1  and  2  when  summed  are  also  equal  to 
the  vector^*.  In  Sector  I,  ^  and^2,  correspond  to  the  states  (P-’n’n)  an  d(p,p,n)  The 
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P  and  n  correspond  to  the  positive  or  negative  IGBT  bus  signals  for  the  respective 
phases  a,  b,  and^7.  The  zero  vectors,  which  are  represented  by  the  states  (P’P’P) 
and(n,n,n '  are  not  shown  but  represented  as  vectors  into  and  out  of  the  page.  The 

magnitude  of  the  vectors  ^  and  2  correspond  to  the  amount  of  time  spent  on  the 

T  T 

switching  states  in  each  sector.  The  duty  cycles  1  and  2  are  the  times  spent  on  each 

V  V  T 

cycle  for  the  vectors  1  and  2 .  The  total  time  spent  for  one  switching  period  is  ' . 


2  TV 
3  T 


(48) 


2  TV 

V.=^£jc_ 

3  T 


(49) 


Applying  the  law  of  sines  to  any  one  of  the  sectors  in  the  SVM  Hexagon  will  produce 


2V*  _  Vj  _  V2 
V3  -sin(6O°-0)-sin(0) 

Substituting  (48  and  49)  into  (50)  is  used  to  find  the  time  of  the  duty  cycles  for  the 
vectors  Vx  and  V2 . 


Tr  = 


'LATssm(m- -e) 

"  dc 


(51) 


T  = 
J  2 


Ts  sin(6») 


dc 


Ts=Ti+T2+T0 


(52) 

(53) 


SVM  has  the  ability  to  minimize  switching  loss  and  harmonic  distortion  by  controlling 
the  order  in  which  the  states  are  applied  [8].  The  switching  pattern  for  each  state  and 

sector  used  in  this  thesis  is  shown  in  Table  2.  The  SVM  algorithm  is  implemented  inside 
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the  SDC  center  where  it  is  oversampled  to  converge  on  the  performance  of  a  true 
continuous  sinusoidal  signal.  The  digital  implementation  of  the  SVM  scheme  used  for 
this  thesis  is  discussed  in  detail  in  [9]  and  the  schematic  diagrams  can  be  found  in 
Appendix  C.  Figure  18  is  a  block  diagram  of  the  functions  for  each  process  in  the  SVM 
block. 
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Table  2.  Space  Vector  Modulation  Switching  Pattern  (From  [9]). 
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Figure  18.  SVM  Digital  Implementation  (From  [9]). 

In  Conclusion,  this  chapter  described  the  methodology  used  to  control  the  DFIG 
for  simulation  and  hardware  design.  The  next  chapter  discusses  the  test  performed  to 
verify  the  operation  of  the  system. 
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IV.  RESULTS  AND  ANALYSIS 


A.  VOLTAGE  AND  CURRENT  COMPARISON 

Several  tests  were  performed  to  verify  the  operation  of  the  wind  energy 
conversion  system.  The  first  test  is  used  to  verify  that  the  system  can  produce  rated  stator 
voltage  and  frequency  on  the  grid.  Figure  19  is  a  comparison  of  the  physical  system  to 
the  simulation  with  the  rotor  current  set  at  1A.  The  simulation  closely  matches  the 
physical  system  and  produces  a  170V  peak  to  peak  signal  at  60  Hz.  The  variables  Vs, 
and  Vs,sim  represent  the  stator  voltage  and  stator  simulation  voltage.  The  variable  Ir 
represents  the  rotor  current.  The  variable  Pout  represents  the  average  instantaneous 
output  power  from  one  phase  multiplied  by  3  delivered  from  the  stator  tenninals  to  the 
grid.  This  notation  is  followed  for  all  the  graphs  in  the  results  section. 


Comparison  of  Stator  Voltages 
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Figure  19.  Comparison  of  Grid  Stator  Voltage  to  Simulation  Source. 
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The  second  test  was  to  verify  the  operation  of  the  system  for  various  changes  in 
rotor  current.  The  control  algorithm  allows  the  user  to  specify  the  amplitude  of  the  rotor 
current.  The  benefits  of  controlling  the  amplitude  of  the  current  are  discussed  in  the  next 
section.  This  test  is  just  to  prove  that  the  system  functions  similarly  to  the  simulation  for 
changes  in  rotor  current.  The  results  are  shown  in  Figure  20. 


Generator  Stator  Currents 
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Figure  20.  Comparison  of  Stator  Currents  to  the  Stator  Simulation  Currents. 

From  Figure  20,  one  can  see  that  as  the  rotor  current  is  increased  in  amplitude,  the  stator 
current  reduces  in  amplitude.  The  simulation  does  not  produce  as  much  change  in  stator 
current  amplitude  as  the  physical  system  for  equal  changes  in  rotor  current.  This  is  not  a 
problem;  the  system  still  produces  voltage  and  current  at  60Hz. 
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B. 


POWER  FACTOR  AND  OUTPUT  POWER  COMPARISONS 


The  next  test  demonstrates  the  usefulness  of  having  the  ability  to  control  the 
amplitude  of  the  rotor  currents.  A  Wind  Energy  Conversion  System  can  only  output  as 
much  power  as  the  wind  inputs  into  to  the  turbine  blades.  A  specific  amount  of  torque  on 
the  turbine  blades  will  generate  a  certain  amount  of  apparent  power  in  the  machine. 
Controlling  the  amplitude  of  the  rotor  currents  allows  the  user  to  detennine  how  much 
real  power  is  abstracted  from  the  total  amount  of  apparent  power  generated  in  the 
machine.  This  technique  is  known  as  power  factor  control  and  shown  in  Figure  21.  The 
voltage  and  current  values  below  are  normalized  to  the  DFIG  ratings. 


Power  Factor  Comparison 


CD 

T3 

=3 

■q. 

E 

< 


CD 

T3 

=3 

■q. 

E 

< 


0.02  0.04 

Time  (s) 

Power  Factor  Comparison 


Vs(lr=3A) 

ls(lr=3A) 


7A 


1 

l  12_ i 

0.02  0.04 

Time  (s) 


CD  . 

-o  1 

Q_ 

<  o 


-i 


Power  Factor  Simulation  Comparison 

2  |- - 1  1 

Vs(lr=1  A) 
ls(lr=1  A) 


CD  . 

tj  1 

Q_ 

<  o 


0  0.02  0.04 

Time  (s) 

Power  Factor  Simulation  Comparison 

2h - 11 

-  Vs(lr=3A) 
ls(lr=3A) 


-1  L 


0.02  0.04 

Time  (s) 


Figure  2 1 .  Power  Factor  Comparison. 

Figure  21  compares  the  power  factor  adjustment  for  the  physical  system  with  the 
simulation.  The  rotor  current  is  adjusted  from  1  to  3  amps  and  the  stator  voltage  and 
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currents  are  graphed.  One  can  see  that  the  angle  between  the  voltage  and  current  are 
almost  in  phase,  as  the  rotor  current  is  shifted  from  1  to  3  amps. 

Figure  22  is  a  display  of  how  the  power  factor  affects  real  power  delivered  to  the 
grid  for  changes  in  rotor  current.  The  power  factor  for  rotor  currents  of  1A,  2A,  and  3 A 
respectively  are  calculated  as  0.3,  0.71,  and  0  .94  for  the  physical  system  and  0.28,  0.69, 
and  0.97  for  the  simulation. 
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Figure  22.  Comparison  of  Output  Power  for  Changes  in  Rotor  Current 

The  last  test  is  used  to  verify  that  the  system  can  operate  under  various  changes  in 
load.  A  DC  machine  is  used  to  simulate  a  change  in  wind  speed  by  increasing  the  amount 
of  torque  delivered  to  the  DFIG  shaft.  A  step  change  of  1  Newton  is  applied  by  the  DC 
drive,  which  is  equivalent  to  175  watts  at  rated  speed.  The  results,  shown  in  Figure  23, 
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prove  that  the  system  is  capable  of  delivering  constant  power  to  the  grid  as  wind  speed 
changes.  The  physical  system  produces  slightly  less  output  power  than  the  simulation  due 
to  mechanical  coupling  losses  from  the  DC  drive  to  the  DFIG. 


Step  Change  in  Torque 


Time  (s) 

Step  Change  in  Torque  Simulation 
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Figure  23.  Comparison  for  Step  Change  in  Torque 

This  chapter  covered  the  results  for  each  test  performed  that  compared  the 
operation  of  the  simulation  to  physical  system.  The  next  chapter  will  discuss  the  results 
and  conclusions. 
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V.  CONCLUSIONS  AND  SUGGESTIONS 


A.  CONCLUSIONS 

A  Wind  Energy  Conversion  System  was  designed,  simulated,  and  constructed 
using  the  methods  discussed.  A  series  of  tests  were  conducted  to  validate  the  design.  The 
system  was  designed  to  work  under  various  load  conditions,  which  in  this  scenario 
simulates  changes  in  wind  speed.  The  ability  to  control  the  rotor  currents  provides  the 
user  with  functionality  to  adjust  power  factor  thereby  increasing  the  efficiency  of  power 
delivered  to  the  grid.  The  response  of  the  Wind  Energy  Conversion  system  was  compared 
to  the  simulation.  The  results  show  that  the  physical  system  behaved  as  predicted  by  the 
simulation. 

A  successful  Wind  Energy  Conversion  system  demonstrates  the  ability  to 
transform  mechanical  energy  delivered  by  the  wind  into  electrical  energy  that  can  be  used 
to  power  any  electrical  grid.  The  Wind  Energy  Conversion  System  successfully 
transformed  mechanical  torque  on  the  shaft  into  electrical  power. 

B.  FUTURE  RESEARCH  OBJECTIVES 

This  thesis  dealt  with  simulation  and  operation  of  a  Wind  Energy  Conversion 
System.  Future  research  will  be  on  various  control  strategies  to  improve  overall 
efficiency. 

One  area  of  research  currently  being  explored  is  the  ability  to  extract  power  from 
both  the  rotor  and  stator  windings  when  the  machine  is  operating  at  super-synchronous 
speed.  This  design  will  require  two  VSI  and  control  software  that  will  allow  for  bi¬ 
directional  flow  of  current  through  the  rotor  windings  for  sub-synchronous  and  super- 
synchronous  operation.  Delivering  power  from  the  rotor  and  stator  windings  will  improve 
the  overall  efficiency  of  the  system. 

The  other  area  of  research  currently  being  explored  is  a  continuation  of 
controlling  the  power  factor  angle.  The  rotor  side  converter  will  be  used  to  control  active 
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and  reactive  power  of  the  DFIG  by  controlling  the  amplitude  of  the  rotor  currents.  The 
amplitude  of  the  rotor  currents  will  follow  a  tracking  characteristic  that  adjusts  generator 
speed  for  optimal  power  generation  depending  on  wind  speed. 
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APPENDIX  A:  DATASHEETS 


SEMISTACK  -  IGBT 


Three-phase  rectifier  + 
inverter  with  brake 
chopper 

SEMITEACH  -  IGBT 
SKM  50  GB  123D 
SKD51 
P3/250F 

Features 


Circuit 

Irms  (A) 

Vgc  /  Vocmax 

Types 

B6CI 

30 

440  /  750 

SEMITEACH  -  IGBT 

Symbol 

Conditions 

Values 

Units 

t-  : 

no  overtoad 

30 

A 

IGBT  -  4x  SKM  50  GB  t23D 

VcES 

1200 

V 

VcEISAT) 

Vg£S 

lc=  50A,  Vge=  15V.  chp  level:  Tj=  25(125)°C 

2,7  (3,5) 
*20 

V 

V 

Ic 

Lcm 

To«»=25  (80)°C 

T„«=25  (90)°C:  1ms 

50(40) 
tOO  (SO) 

A 

A 

Hectitier  -  lx  SkD  5f/14 

without  filter 

3x480 

V 

with  filter 

3x380 

V 

VDOn« 

DC  Capacitor  bank  -  Electrolytic  2x  2200pF/400V 

total  equivalent  capacitance 

max.  DC  voltage  applied  to  the  capacitor  bank 

1100/800 

750 

pF/V 

V 

Driver  -  4x  SKHI  22 

Power 

supply 

Current 

0/15 

V 

con  sump 
tion 

max:  per  driver 

16 

mA 

Thermal 

trip 

Normally  Open  type  (NO) 

71 

°C 

•  Multi-function  IGBT  converter 

•  Transparent  enclosure  to  allow 
visualization  of  every  part 

•  IP2x  protection  to  minimize 
safety  hazards 

•  External  banana/BNC  type 
connectors  for  all  devices 

•  Integrated  drive  unit  offerng 
short-circuit  detection/cut-off, 
power  supply  failure  detection, 
interlock  of  IGBTs  +  galvanic 
isolation  of  the  user 

•  Forced-air  cooled  heatsink 


Typical  Applications 

•  Education:  One  stack  can 
simulate  almost  all  existing 
industrial  applications: 

-  3-phase  inverter+brake  chopper 

-  Buck  or  boost  converter 

-  Single  phase  nverter 

-  Single  or  3-phase  rectifier 


'I  Photo  non-contract Ljal 


This  technical  information  specifies  semiconductor  devices  but  promises  no 
characteristics.  No  warranty  or  guarantee  expressed  or  implied  is  made  regarding 
delivery,  performance  or  suitability. 
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Model  8231  Three-Phase  Wound-Rotor  Induction 


Each  phase  of  the  stator  windings  of  this  motor  is  inde¬ 
pendently  terminated  and  identified  on  the  faceplate  to 
permit  operation  in  either  delta  or  star  (wye)  configura¬ 
tion.  The  rotor  windings  are  brought  out  to  the  faceplate 
via  external  slip  rings  and  brushes.  This  machine  can  be 
used  as  a  wound-rotor  induction  motor,  phase  shifter, 
single-phase  variable  coupling  transformer,  three-phase 
transformer,  selsyn  control,  frequency  converter  or 
asynchronous  induction  generator.  The  speed  of  this 
machine  can  be  controlled  through  the  use  of  the  Three- 
Phase  Rheostat  (Model  8731). 


SPECIFICATIONS 


Model  821 1  DC  Motor/Generator 

120/208  V-  60  Hz 

220/380  V-  50  Hz 

240/415  V- 50  Hz 

Power  Requirement 

120/208  V 

220.380  V 

240415  V 

Rating  Motor  Output  Power 

Generator  Output  Power 

Armature  Voltage 
Shunt  Field  Voltage 

Full  Load  Speed 

Ful  Load  Motor  Current 

Full  Load  Generator  Current 

175  W 

120  W 

now 

120  W 

120  V -DC 

220  V-  DC 

240  V- DC 

120  V-  DC 

220  V-  DC 

240  V- DC 

ISOOrimin 

1500  r/mn 

1500  r/mn 

2.8  A 

ISA 

1.1  A 

1  A 

0.5  A 

0.5  A 

Physcal  Characteristics  Dimensions  (H  x  W  x  D) 

Net  Weight 

308  x  291  x  440  rmi  (12.1  x  11.5  x  17.3  n) 

14.1  kg  (31  b) 

Model  8221  Four-Pole  Squirrel-Cage  Induction  Motor 

120/208  V-  60  Hz 

220/380  V- 50  Hz 

240/415  V-  50  Hz 

Power  Requirement 

120/208  V 

220.380  V 

240415  V 

Ratng  Output  Power 

Stator  Voltage 

Full  Load  Speed 

Full  Load  Current 

175  W 

120/208  V,  3-phase 

220/360  V.  3-phase 

240415  V,  3-phase 

1670  r/mm 

1360  r/mn 

1395  r/mn 

1.2  A 

0.52  A 

0  46  A 

Physical  Characteristics  Dimensions  (H  x  W  x  D) 

Net  Weight 

308  x  201  x  440  mm  ( 12.1  x  1 1.5  x  17.3  in) 

13.5  kg  (29  7  lb) 

Model  8231  Three-Phase  Wound-Rotor  Induction  Motor 

120/208  V-  60  Hz 

220/380  V  -  50  Hz 

240/415  V- 50  Hz 

Power  Requirement 

120/208  V 

220.380  V 

240415  V 

Ratng  Output  Power 

Stator  Voltage 
Rotor  Voltage 

Full  Load  Speed 

Full  Load  Current 

175  W 

120/208  V,  3-phase 

220/380  V.  3-phase 

240415  V,  3-phase 

60/1 04  V,  3-phase 

110/190  V.  3-phase 

1201208  V.  3-phase 

1500nVnin 

1240  r/mn 

1315  r/mn 

1.3  A 

0.53  A 

0.48  A 

Phys.cal  Characteristics  Dimensions  (H  x  W  x  D) 

Net  Weight 

308  x  291  x440  mm  (12.1  x  11.5x  17.3  in) 

14  kg  (30.8  lb) 

Model  821 1  DC  Motor  Generator 


This  machine  can  be  run  independently  as  a  DC  motor 
or  a  DC  generator.  The  armature,  shunt  field,  and  series 
field  windings  are  terminated  separately  on  the  faceplate 
to  permit  long  and  short  shunt  as  well  as  cumulatively 
and  differentially  compounded  motor  and  generator  con¬ 
nections.  This  machine  is  fitted  with  exposed  movable 
brushes  to  allow  students  to  study  the  effect  of  armature 
reaction  and  commutation  whi  le  the  machine  is  operating 
under  load.  An  independent  circuit-  breaker  protected, 
shint-field  rheostat  is  mounted  on  the  faceplate  for  motor 
speed  control  or  generator  output  voltage  adjustment. 


POWER  SUPPLY 
MODEL  8821 


The  Power  Supply  provides  fixed  and  variable  AC 
and  DC  voltage  sources,  all  terminated  by  color-coded 
4  mm  safety  sockets.  Independent  circuit  breakers,  reset 
at  the  front  panel,  protect  the  input  to  and  output  from  the 
Power  Supply.  Indicator  lamps  monitor  the  presence  of 
input  voltage  in  each  phase.  When  a  phase  leg  of  the 
site’s  power  service  is  out,  the  lamp  goes  off  to  reflect  this 
condition. 


A  voltmeter,  connected  through  a  selector  switch, 
monitors  the  variable  AC  and  DC  outputs  and  fixed  DC 
output  A  24  V  AC  output  provides  a  low-voltage  supply 
required  to  operate  other  EMS  equipment  such  as  meter¬ 
ing  modules  and  modules  used  in  the  Power  Electronics 
Training  System. 


SPECIFICATIONS 


Model  8821  -  Power  Supply 

1201208  V-  60  Hz 

220/380  V  -  50  Hz 

240/415  V  -  50  Hz 

Input  Line  Voltage 

Line  Current 

Service  Installation 

120i’208V 

220/380  V 

240/415  V 

15  A 

10  A 

20  A,  3-phase.  5  wres.  star  (Wye)-  connected,  including  neutral  and 
ground 

Outputs  Fixed  AC  3-Phase 

Variable  AC  3-Phase 

Vanable  DC 

Fated  DC 

Low  Power  AC 

120/208  V-  15  A 

220/380  V-  10  A 

240/415  V-  10  A 

0-120/208  V- 5  A 

0-220/380  V-3  A 

0-240/415  V-3  A 

0-120  V-8  A 

0-220  V  -  5  A 

0-240  V  —  5  A 

120V-2A 

220  V  -  1  A 

240  V- 1  A 

24  V  -  3  A 

Wall  Outlet  included 

NEMA  L21-20 

NEMA  L22-20 

CLIPSAL  5630520 

Power  Cord 

3  m  (10  ft) 

Physical  Characteristics  Dimensions  (HxWx  D) 

Net  Weight 

308  x  287  x  500  mm  (12.1  x  1 1.3  x  19.7  in) 

18  4  kg  (40.5  lb) 
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APPENDIX  B:  MATLAB  M-FILES 


A.  MATLAB  INITIAL  CONDITIONS  FILE 


Dbl  Fed  Ind.m 


clc, clear 
Vdc=150; 

DC_Bus_Prot=200; 
omega_b  =  2*pi*60; 
oversample=2 ; 
pulsect  =  1800/oversample; 
step_ct=l ; 
clock_freq=25e6; 

tstep=.0004;  %  only  to  speed  up  simulation. 

%tstep=step  ct/clock  freq;  %When  compiling  software  us  this  line 


%Bus  voltage  maintained  for  Capacitor  on  VSI 
%Secure  firing  circuit  at  200Vdc. 

%  Base  Frequency  in  radians/second 
%  Determine  oversampling  rate  in  SVM  block 


%More  constants  used  for  simulation 
twopiby3  =  2*pi/3; 
poles  =  4; 

polesby2 J=poles* 12 /2 / .089; 

TL= . 25*746/ (2/poles*omega_b) ; 


%Calculated  constants  from  Induction  motor  test. 
rs=12 ; 
rr  =  15; 

XI s  =  9.1; 

Xm  =126; 

Xlr  =9.1; 


%D= (Xls+Xm) * (Xlr+Xm) -XmA2 ;  Defined  in  book  for  per  unit  values  of 
Inertia . 

rsbyXls  =  rs/Xls; 
rrbyXlr  =  rr/Xlr; 

Xaq  =  1/ ( 1/Xm+l/Xls+l/Xlr) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

XaqbyXm  =  Xaq/Xm; 

XadbyXm  =  Xad/Xm; 

psi_qsic=0 ; 
psi_dsic=0 ; 
psi_qric=0 ; 
psi_dric=0 ; 

omegar  ic  =omega_b;  %Set  the  initial  condition  speed  of  rotor  to  1800 
rpm. 

Kp  i=.231;  %current  Proportional  gain 
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Ki  i=415 . 7*1 . 042;  %Current  Integral  gain 

%Used  for  Mealy  State  Machine  in  A/D  Current  simulation. 

F_mat  =  [0  0  0  1 ; 1  1  2  0;2  2  3  0;3  3  0  0 ]  ; 

0  mat  =  F  mat; 


%Everything  below  this  line  is  used  for  the  encoder  controls 
%fvtool (outputl , output2 ) 

t  square=2e-4 / 1 0  ;  %2e-4  produces  1500RPM.  . . 4e-4  produces  750RPM 
%8e-4  =  374RPM. . . 10e-4=300RPM 

%there  is  400  1/2  pulses  per  rotation  therefore: 

%60sec/ (400*  l/2*t_square) =  correct  RPM 
RPM=60/ (200*t_square) ; 

Hz=RPM/60; 

sw_freq=15000; 

%  sw  counter=round ( f  clock/sw  f  req-mod  ( f__clock/sw  freq, 10)); 

%Counter  for  sawtooth  for  switching  modulo  10  used  so  step_ct  can  be  10 
Total_Rotations=6; 

%Ctr=  [  1 :  2  A8 ]  /2''8/tstep/4  00*50;  %16,  000,000 
%reciprocal=l ./Ctr;%l/16, 000, 000=2 44. 1406 
Ctr= [ 1 : 2  A11 ]  ; 

reciprocal=3 60/30/1 ./Ctr;%360/30/l./Ctr 
%  [Pk,  Nk,  0k] 

output_vec= [ 0 ; . . . %Nothing  is  occurring  with  [0,0,0] 

5;...%3Motor  is  turning  in  the  CCW  direction  and  crossing 

the 

%zero  pt  with  [0,1,1] 

0 ; . . . %0Nothing  is  occurring  with  [0,0,0] 

4;...%2Motor  is  turning  in  the  CCW  direction  with  [0,1,0] 

0 ;  .  .  .  %0Nothing  is  occurring  with  [0,0,0] 

3;  .  .  . %5  [1, 0, 1]  indicates  CW  rotation  and  resets  rotation 
%counter 

0 ; . . . %0Nothing  is  occurring  with  [0,0,0] 

2  ;  .  .  .  %4  [  1 , 0 , 0 ]  indicates  CW  rotation 

3; . . . %5 [1, 0, 1]  indicates  CW  rotation  and  resets  rotation 
%counter 

0 ;  .  .  . %0Nothing  is  occurring  with  [0,0,0] 

2 ; . . . %4 [ 1 , 0 , 0 ]  indicates  CW  rotation 
0 ;  .  . . %0Nothing  is  occurring  with  [0,0,0] 

5;  . . . %3 [0, 1, 1]  indicates  CCW  direction  and  resets  rotation 
%counter 

0 ; . . . %0Nothing  is  occurring  with  [0,0,0] 

4;...%2  Motor  is  turning  in  the  CCW  direction  with  [0,1,0] 
0];  %0  Nothing  is  occurring  with  [0,0,0] 


next 

output=next; 


[0, 0,1, 0;1, 0,1,0]; 
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B 


MATLAB  M-FILE  USED  FOR  SPACE  VECTOR  MODULATION 


cverflow3.m 

function  [sectorl,  sector2,  sector3,  sector4,  sector5, 
overflow3 (x) 

%gain  =  xfix ( {xlUnsigned, 10, 7 } , 2 . 359296/3) ; %for  60  hz 

gain  =  xfix ( {xlUnsigned,  10, 7 },  2 . 359296)  ;  %for  180  hz 

%tempv=gain*x; 

tempv=x ; 

if  tempv<=171-l 

sectorl=xfix ( {xlBoolean} ,  1)  ; 
sector2=xfix ( {xlBoolean} ,  0) ; 
sector3=xfix ( {xlBoolean} , 0) ; 
sector4=xfix ( {xlBoolean} , 0) ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xfix ( {xlUnsigned, 10,0}, tempv) ; 
elseif  tempv<=2*171-l 

sectorl=xfix ( {xlBoolean} ,  0) ; 
sector2=xf ix ( {xlBoolean } , 1 ) ; 
sector3=xfix ( {xlBoolean} , 0) ; 
sector4=xfix ( {xlBoolean} , 0) ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xf ix ( { xlUnsigned, 10,0}, tempv- 17 1 ) ; 
elseif  tempv<=3*171-l 

sectorl=xfix ( {xlBoolean} , 0) ; 
sector2=xfix ( {xlBoolean} ,0); 
sector3=xfix({xlBoolean},l) ; 
sector4=xfix ( {xlBoolean} , 0) ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xf ix ( { xlUnsigned, 10,0}, tempv- 2  *  171 )  ; 
elseif  tempv<=4*171-l 

sectorl=xfix ( {xlBoolean} ,  0)  ; 
sector2=xfix ( {xlBoolean} , 0) ; 
sector3=xfix ( {xlBoolean}  ,0); 
sector4=xfix ({ xlBoolean }, 1 ) ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xfix ( {xlUnsigned, 10,0}, tempv-3*171) ; 
elseif  tempv<=5*171-l 

sectorl=xfix ( {xlBoolean} , 0) ; 
sector2=xfix ( {xlBoolean} , 0) ; 
sector3=xfix ( {xlBoolean} , 0) ; 
sector4=xfix ( {xlBoolean} ,0); 
sector5=xfix({xlBoolean},l) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xf ix ( { xlUnsigned, 10,0}, tempv- 4  *  17 1 )  ; 

else 


sector6. 


sectorl=xfix ( {xlBoolean} , 0) ; 
sector2=xfix ( {xlBoolean} ,  0)  ; 
sector3=xfix ( {xlBoolean} , 0) ; 
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end 


sector4=xfix ( {xlBoolean} , 0) ; 
sector5=xfix ( {xlBoolean} ,  0)  ; 
sector6=xf ix ({ xlBoolean },  1 ) ; 
z=xf ix ( { xlUnsigned, 10,0}, tempv-5* 17 1 )  ; 


ramp2mod.m 

function  z  =  ramp2 (x) 

gain=xfix ( {xl Signed, 20, 19}, 1/1800) 

z=xf ix ( { xl Signed,  14  ,  13 } ,  x*gain)  ; 


thetaconv2.m 

function  [y]  =  thetaconv(x) 
gainl  =  xf ix ({ xlSigned, 14 , 1 0 },  2*3 . 14  )  ; 
gain2  =  xf ix ({ xlSigned, 14 , 1 0 }, 1 /gainl ) 
if  x<0 

y=xf ix ( { xlUnsigned, 10,0}, (x+gainl ) *gain2* 1024 ) ; 
else 

y=xf ix ( { xlUnsigned, 10,0}, x*gain2*1024  )  ; 
end 


C.  MATLAB  M-FILES  USED  FOR  ENCODER 
mcode5A.m 

function  [NewRotation,  NewValue,  Direction]  =  mcode5A(Pk,  Nk,  Zk, 
OldRotation,  OldValue) 

%This  fuction  records  CW  and  CCW  rotations  via  matlab  MCODE  simulation 
%block.  It  also  records  360  degrees  of  rotation  and  total  degrees 
traveled 

%Author:  LT  Andrew  M  LaValley 
%Last  Modified:  14  APR  08 

if  Pk  &&  ~Zk  &&  ~Nk  %Increments  the  #  of  pulses  produced  in  the  CW 
rotation 


Direction=0 ; %indicates  positive  direction 

if  01dValue==399  %this  keeps  the  value  positive  if  the  value  is 
zero . 


NewValue=xf ix ( {xlSigned,  15, 0 } ,  0)  ; 

NewRotation=xf ix ( { xlSigned, 15,0},  OldRotation+1 )  ; 


NewValue=xf ix ( { xlSigned, 15,0},  OldValue+1 )  ; 
NewRotation=xf ix ( { xlSigned, 15,0}, OldRotation) ; 
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end 

elseif  Nk  &&  ~Zk  &&  ~Pk  %Decrements  the  #  of  pulses  produced  in  the  CCW 
rotation  if  not  zero 

Direction=l ; %Indicates  negitive  direction 

if  OldValue==0  %this  keeps  the  value  positive  if  the  value  is  zero. 

NewValue=xf ix ( {xl Signed,  15, 0 },  399)  ; 

NewRotation=xf ix ( { xl Signed, 15,0},  OldRotation-1 )  ; 

%Direction=l ; 

else 

NewValue=xf ix ( { xl Signed, 15,0}, OldValue-1 ) ; 

NewRotation=xf ix ( { xl Signed, 15,0}, OldRotation) ; 

end 


else  %01dValue  will  equal  NewValue  and  OldRotation  will  equal  NewRotion 
if 

%the  above  if/else  statements  do  not  apply. 

Direction=0 ; %def ault  direction 
NewValue=xf ix ( {xlSigned, 15,0}, OldValue) ; 

NewRotation=xf ix ( { xlSigned, 15,0}, OldRotation) ; 

end 


D.  MATLAB  M-FILES  FOR  CHIPSCOPE  INTERFACE 
Black_box_dc_machine2_config.m 

function  code  config(this  block) 

%  Revision  History: 

o, 

o 

%  18-Dec-2008  (15:15  hours): 

%  Original  code  was  machine  generated  by  Xilinx's  System 

Generator 

%  after  parsing  H:\Docs\work  f iles\Docs\classes\EC3130  2009\DC 

%  machine  lab\black  box  dc  machine2.vhd 

this  block. setTopLevelLanguage ( ' VHDL ' ) ; 

this  block. setEntityName ( ' code ' ) ; 

%  System  Generator  has  to  assume  that  your  entity  has  a 
combinational 
%  feed  through; 

%  if  it  doesn't,  then  comment  out  the  following  line: 

45 


this  block . tagAsCombinational ; 


this  block. addSimulinklnport ( ' ind ' ) ; 
this  block . addSimulinklnport (' ila  clock'); 
this  block. addSimulinklnport ( ' ind2 ' ) ; 

this  block. addSimulinkOutport ( ' outd ' )  ; 
this  block . addSimulinkOutport (' open  loop'); 
this  block. addSimulinkOutport ( ' speed ' ) ; 

outd_port  =  this_block . port (' outd ') ; 
outd_port . setType ( ' UFix_l_0 ' ) ; 

open_loop_port  =  this_block . port ( ' open_loop ' ) ; 
open_loop_port . setType ( ' UFix_l_0 ' ) ; 
speed_port  =  this_block . port (' speed ') ; 
speed_port . setType ( ' UFix_8_0 ' ) ; 


o,  _ 

o 

if  (this  block . inputTypesKnown) 

%  do  input  type  checking,  dynamic  output  type  and 
%  this  code  block. 

if  (this  block . port (' ind '). width  ~=  1); 

this  block .  setError  (' Input  data  type  for  port 
width=l . ' ) ; 
end 

this  block. port ( ' ind' ) . useHDLVector (false) ; 

if  (this  block . port (' ila  clock '). width  ~=  1); 

this  block .  setError  (' Input  data  type  for  port 
have  width=l . ' ) ; 
end 

this  block . port (' ila  clock' ) .useHDLVector (false) ; 

if  (this  block . port (' ind2 '). width  ~=  48); 

this  block .  setError  (' Input  data  type  for  port 
width=4  8 . ' ) ; 
end 

end  %  if ( inputTypesKnown) 

O,  _ 

o 


if  (this  block . inputRatesKnown) 

setup  as  single  rate (this  block, 'elk' 
end  %  if ( inputRatesKnown) 


'  ce ' 


%  Add  addtional  source  files  as  needed. 
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generic  setup  in 


"ind"  must  have 


"ila  clock"  must 


"ind2"  must  have 


%  | -  |  Add  files  in  the  order  in  which  they  should  be 

%  compiled.  |  If  two  files  "a.vhd"  and  "b.vhd"  contain  the  entities 
%  entity  a  and  entity  b,  and  entity  a  contains  a  |  component  of  type 
%  entity_b,  the  correct  sequence  of  |  addFileO  calls  would  be: 

%  this  block . addFile (' b . vhd ') ;  |  this  block . addFile ( ' a . vhd ' ) ; 


%  this  block . addFile ('') ;  this  block . addFile ('') ; 
this  block . addFile (' black  box  dc  machine2 . vhd ' ) ; 

return ; 


function  setup  as  single  rate (block, clkname, cename) 
inputRates  =  block . inputRates ; 
uniquelnputRates  =  unique ( inputRates ) ; 

if  ( length (uniquelnputRates ) ==1  &  uniquelnputRates ( 1 ) ==Inf) 

block . setError (' The  inputs  to  this  block  cannot  all  be  constant. 1 ) ; 
return; 
end 

if  (uniquelnputRates (end)  ==  Inf) 
hasConstantlnput  =  true; 

uniquelnputRates  =  uniquelnputRates ( 1 : end-1 ) ; 

end 

if  ( length (uniquelnputRates )  ~=  1) 

block . setError (' The  inputs  to  this  block  must  run  at  a  single 
rate .  '  )  ; 

return; 

end 

thelnputRate  =  uniquelnputRates ( 1 ) ; 
for  i  =  1 : block . numSimulinkOutports 

block. outport (i) . setRate (thelnputRate) ; 

end 

block . addClkCEPair (clkname, cename, thelnputRate) ; 
return; 
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APPENDIX  C:  SIMULINK/  XILINX  MODEL  OF  WIND  ENERGY 

CONVERSION  SYSTEM 
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Dbl_Fed  _lnd_6/A  toD  conversion  (current) 
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Dbl  Fed  lnd_6/A  toD  conversion  (current)ZSubsystem 
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Dbl_Fed_lnd_6/A  toD  conversion  (current)/Subsystem/pulse1 
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Dbl_Fed  _lnd_6/A  toD  conversion  (voltage) 
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Dbl_Fed_lnd_6/Chipscope  interface 
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Dbl  Fed  Ind  6/Controller 
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APPENDIX  D:  TRANSFORMATION  DERIVATION 


vab  ~va~vb’vbc~vb~vc~vb~(~va~vb)<  where  va  +  v6  +  =  0  for  a  wye  connected  capacitor  bank  (with  the  neutral  floating ) 

Note  also  that  \a  +  ib  +4  =  0  for  a  wye  connected  transformer  winding  set. 
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